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INTRODUCTION 

Spinal  cord  injuries  (SCI)  produce  direct  mechanical  disruption  with  subsequent  severe 
degeneration  of  axons,  and  are  the  processes  underlying  the  associated  neurologic  deficits 
observed  in  such  injuries.  Histological  studies  of  fixed  tissue  in  animal  models  of  SCI  have 
described  axonal  loss  and  demyelination  occurring  after  trauma.  Research  at  the  University  of 
Pennsylvania  site  brings  novel  magnetic  resonance  methodology  to  bear  with  the  objective  of 
obtaining  quantitative  information  on  axonal  degeneration  and  myelin  loss  following  spinal  cord 
injury  in  a  mouse  model  by  pursuing  the  following  Specific  Aims  per  the  work  statement: 

1.  We  will  perform  q-space  MR  imaging  (QSI)  and  simulations  of  QSI  to  quantify  axonal 
architecture  in  healthy  and  injured  mouse  spinal  cords. 

2.  We  will  quantify  myelin  content  with  three  quantitative  MRI  techniques  in  healthy  and 
injured  mouse  spinal  cords  and  compare  the  results  with  histology. 

Specific  Aim  1 : 

Subsequent  to  the  upgrade  of  the  Bruker  NMR/MRI  system,  hardware  modifications  were 
needed  to  connect  our  previous  custom  gradient  coil  to  the  new  system.  The  gradient  coil  also 
had  to  be  optimized  and  recalibrated  for  the  new  system.  This  included  a  systematic  and 
thorough  analysis  of  RF  noise  sources  in  the  hardware,  resulting  in  the  implementation  of  new 
shielding  and  a  new  filter  box  on  the  gradient  leads  to  isolate  the  receive  chain  from  both  room- 
ambient  RF  signals  and  RF  noise  leakage  from  the  gradient  amplifier.  Our  implementations 
were  similar  in  design  to  those  used  by  Bruker  for  their  gradient  coils.  In  addition,  we  discovered 
a  crack  in  the  epoxy  of  our  custom  gradient  coil  that  produced  a  vibrational  phase  instability  in 
the  data,  and  we  corrected  this  by  applying  a  new  layer  of  epoxy  resin  to  the  gradient  coil 
windings.  After  these  modifications  to  the  gradient  coil  hardware  had  been  accomplished,  we 
performed  a  new  calibration  of  our  custom  gradient  coil  as  this  was  needed  to  accurately 
calculate  axon  diameters  from  the  q-space  data.  To  ensure  a  reliable  calibration  we  used  three 
independent  methods  that  agreed  with  each  other:  1)  a  high-resolution  transverse  image  of  a 
capillary  tube  of  known  inner  diameter  (measured  by  micro-CT)  filled  with  PEG  in  water  was 
acquired  using  our  custom  gradient  for  the  frequency  axis.  The  Bruker  gradient  for  the  phase 
axis  and  the  frequency-axis  scale  factor  was  adjusted  to  give  equal  frequency/phase  tube 
diameters;  2)  projections  of  the  same  PEG/water  tube  were  acquired  orthogonal  to  the  tube  at 
varying  gradient  amplitudes  using  our  custom  coil  and  the  resulting  NMR  line  widths  were 
plotted  versus  actual  current  applied  to  the  coil;  and  3)  the  ADC  of  deionized  water  was 
measured  using  our  custom  coil  and  compared  to  the  literature  value  at  the  same  temperature. 
Finally,  as  the  old  QSI  pulse  sequence  program  did  not  run  on  the  new  system,  a  new  QSI 
pulse  sequence  program  had  to  be  developed. 

Once  the  new  pulse  sequence  program  has  been  tested,  QSI  experiments  will  be 
performed,  followed  by  histologic  analysis  and  QSI  simulations. 

Progress  has  been  made  toward  translation  of  the  QSI  methodology  to  the  clinic  (as 
reported  in  part  previously).  The  pulse  sequence  designed  toward  the  end  of  the  prior  reporting 
year  has  since  yielded  QSI  displacement  maps  on  1.5T  on  a  clinical  imager  (Siemens  Sonata) 
on  fixed  pig  spinal  cords.  Since  the  pig  spinal  cord  is  similar  in  size  to  the  human  spinal  cord 
and  the  experiments  were  performed  with  standard  imaging  gradients  the  results  suggest 
feasibility  of  performing  studies  in  humans  with  spinal  cord  injury  (even  though  this  was  not  a 
specific  objective  of  the  current  project).  This  research  was  presented  as  a  poster  at  the  2012 
Annual  Meeting  of  the  International  Society  for  Magnetic  Resonance  in  Medicine  in  Melbourne, 
Australia.  Notably,  we  have  been  able  to  distinguish  white  matter  tracts  differing  in  axon 
diameter  and  density.  Figure  la  shows  analyzed  white  matter  regions.  Figure  1b  shows  the 
mean  displacements  in  the  three  regions,  commensurate  with  the  different  axon  sizes  (1). 
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Figure  1  a)  Sample  transverse  images  of  pig  spinal  cord  acquired  at  1.5T  with  regions  of 
interest  indicated  (white  bar  =  5mm);  b)  RMS  displacement  obtained  from  QSI  displacement 
maps  consistent  with  known  axon  diameters  in  the  three  WM  tracts  examined. 


Specific  Aim  2: 

We  previously  reported  preliminary  results  on  our  endeavors  toward  direct  imaging  of 
myelin  in  the  spinal  cord  of  rodents  by  3D  ultra-short  echo-time  (UTE)  MRI.  This  work  has  now 
been  published  in  the  Proceedings  of  the  National  Academy  of  Sciences  (PNAS)  (2).  In  brief: 
We  identified  the  spectrum  of  myelin  in  the  spinal  cord  in  situ  as  well  as  in  myelin  lipids 
extracted  via  a  sucrose  gradient  method,  and  investigated  its  spectral  properties.  High- 
resolution  solution  NMR  spectroscopy  showed  the  extract  composition  to  be  in  agreement  with 
myelin's  known  chemical  make-up.  The  400-MHz  spectrum  of  the  myelin  extract,  at  20°C 
(room  temperature)  and  37°C,  consists  of  a  narrow  water  resonance  superimposed  on  a  broad 
envelope  shifted  approximately  3.5  ppm  upfield,  suggestive  of  long-chain  methylene  protons. 
Superimposed  on  this  signal  are  narrow  components  resulting  from  functional  groups  matching 
the  chemical  shifts  of  the  constituents  making  up  myelin  lipids.  The  spectrum  could  be  modeled 
as  a  sum  of  super-Lorentzians  with  a  T2*  distribution  covering  a  wide  range  of  values  (0.008-26 
ms).  Overall,  there  was  a  high  degree  of  similarity  between  the  spectral  properties  of  extracted 
myelin  lipids  and  those  found  in  neural  tissue.  The  normalized  difference  spectrum  had  the 
hallmarks  of  membrane  proteins  not  present  in  the  myelin  extract.  Using  3D  radially  ramp- 
sampled  proton  MRI,  with  a  combination  of  adiabatic  inversion  and  echo  subtraction,  the 
feasibility  of  direct  myelin  imaging  in  situ  is  demonstrated.  Last,  the  integrated  signal  from 
myelin  suspensions  is  shown,  both  spectroscopically  and  by  imaging,  to  scale  with 
concentration,  suggesting  the  potential  for  quantitative  determination  of  myelin  density  in  both 
healthy  and  injured  spinal  cords. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Demonstrated  feasibility  of  translation  of  the  Q-Space  mapping  technique  of  axon 
diameter,  intra-  and  extra-axonal  space,  from  the  laboratory  to  the  clinic. 

•  Isolated  myelin  and  showed  that  its  NMR  properties  match  those  of  myelin  in  situ  in  the 
spinal  cord. 

•  Demonstrated  a  path  toward  quantitative  myelin  imaging  in  humans. 

•  Identified  the  magnetic  resonance  signal  of  myelin  and  showed  feasibility  of  direct 
imaging  of  neural  myelin  as  a  new  metric  for  the  evaluation  of  spinal  cord  injury. 

•  Fully  characterized  the  relaxation  properties  of  myelin  toward  a  more  complete 
understanding  of  the  MRI  signal  in  quantitative  myelin  imaging  (Wilhelm  et  al,  20th 
Annual  ISMRM  Scientific  Meeting,  Melbourne,  Australia  (3)) 

•  Published  major  paper  in  the  Proceedings  of  the  National  Academy  of  Sciences 
(Wilhelm,  etal,  PNAS  2012). 
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•  Performed  experiments  on  the  upgraded  Bruker  Instruments  micro-imaging  system  after 
interfacing  custom-built  gradients  for  high-resolution  q-space  imaging  and  achieved 
performance  goals. 

OUTCOMES 

The  new  myelin  imaging  technique  developed  under  this  grant  has  shown  feasibility  for 
quantitative  assessment  of  myelin  content  in  the  CNS  of  the  injured  rat  spinal  cord. 


CONCLUSION 

While  substantial  progress  has  been  made  toward  achievement  of  the  goals,  the  main  objective 
to  obtain  data  in  the  injured  cord  and  the  hypothesized  temporal  changes  has  not  been  attained 
yet,  largely  due  to  the  technical  problems  that  have  now  been  solved.  Our  partners  at  Drexel 
University  have  provided  the  injured,  fixed  spinal  cords  and  we  are  confident  we  will  complete 
the  project  during  year  3  of  this  award. 
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Feasibility  of  low  q-space  diffusion  MRI  at  1.5T 

Henry  H.  Ong1,  Yusuf  Bhagat1,  Jeremy  Magland1,  and  Felix  W.  Wehrli1 

laboratory  for  Structural  NMR  Imaging,  Departement  of  Radiology,  University  of  Pennsylvania  School  of  Medicine,  Philadelphia,  PA,  United  States 

Introduction 

By  exploiting  the  regularity  of  molecular  diffusion  restrictions  such  as  axon  membranes  and  myelin  sheaths1,  q-space  imaging2,3  (QSI)  offers  po¬ 
tential  for  indirect  assessment  of  white  matter  (WM)  axonal  architecture.  For  example,  QSI  can  accurately  estimate  mean  axon  diameter  (MAD)  and 
intracellular  volume  fraction  (ICF)4,5.  Unfortunately,  the  application  of  QSI  on  a  clinical  scanner  is  severely  constrained  by  the  low  gradient  strengths 
available,  which  limits  the  maximum  achievable  q-value  (q  =  (27i)'1yG5,  G  =  gradient  amplitude,  and  5  =  gradient  duration).  Low  maximum  q-value 
leads  to  insufficient  displacement  resolution  to  accurately  study  axons,  which  have  an  MAD  of  1-3  pm.  Low  q-value  diffusion  MRI5,  in  which  axonal 
architecture  information  is  extracted  by  fitting  the  q-space  signal  decay  (E(q))  at  low  q-values  (q'1  »  MAD)  under  the  short  gradient  pulse  approxi¬ 
mation  (SGPA),  does  not  require  high  gradient  amplitudes.  However,  low  clinical  gradient  strengths  lead  to  violation  of  SGPA.  Here,  we  test  the  fea¬ 
sibility  of  implementing  low  q-value  diffusion  MRI  on  a  1.5T  system  by  assessing  axonal  architecture  in  excised  fixed  pig  spinal  cords. 

Methods 

As  described  in  [5],  at  low  q-values  (q  -*»  MAD),  the  signal  decay  is  given  by  E(q)=exp(-27t2q2Z2) 

(Eq.l),  where  Z  is  the  root  mean  squared  (RMS)  displacement  of  diffusing  molecules  during  a  diffusion 
time  A6,7.  As  E(q)  contains  signal  from  extra-  and  intra-cellular  spaces  (ECS  and  ICS),  a  two- 
compartment  version  of  Eq.l  can  be  defined:  E(q)=  fE-exp(-27i2q2ZE2)  +  fiexp(-2n2q2Z2)  (Eq.  2),  where 
fE  and  f\  are  the  relaxation-weighted  ECS  and  ICS  volume  fractions  and  ZE  and  Zx  are  the  RMS  displace¬ 
ment  of  diffusing  molecules  in  the  ECS  and  ICS.  From  Eq.  2,  MAD  and  ICF  can  be  estimated  from  Z\  and 
fi,  respectively. 

For  validation  of  this  method,  five  fixed  cervical  spinal  cords  (SC)  harvested  from  five  skeletally  ma¬ 
ture  Yucatan  mini-pigs  were  used.  Before  experiments,  the  SCs  were  placed  in  tubes  filled  with  Fomblin 
(Sigma-Aldrich)  to  keep  the  specimens  hydrated  and  to  remove  any  background  signal.  The  low  q-value 
diffusion  MRI  method  was  implemented  on  a  1.5T  Siemens  Sonata  MRI  scanner  (Erlangen,  Germany) 
with  40  mT/m  gradients  using  a  custom  single-slice  PGSE  with  multi-shot  fly-back  EPI  readout  pulse  se¬ 
quence.  The  body  coil  was  used  for  transmit  and  a  custom-built  4-channel  phased  array  coil  (Insight  MRI) 
was  used  for  receive.  The  imaging  parameters  were:  A/8/TE  =  98.7/55/257ms,  128x128,  FOV  =  64x64 
mm,  slice  thickness  =10  mm,  number  of  shots  =  8,  NA  =  36,  and  TR  =  2  s.  The  diffusion  gradient  was 
applied  perpendicular  to  the  SCs  in  32  increments  (qmax  =  0.08  pm'1)  and  the  scan  time  -5  hours.  Note 
that  these  values  for  5  and  A  violate  SGPA.  All  five  SCs  were  imaged  simultaneously.  After  Fourier 
transform,  a  3D  matrix  of  32  2D  images  at  various  q-values  was  obtained.  An  average  E(q)  was  measured 
in  ROIs  within  the  dorsal,  ventral,  and  lateral  columns  of  the  SCs  (Fig.  1).  This  average  E(q)  was  fit  with 
Eq.  2  under  the  constraint  fg+f^l. 

Results  and  Discussion 

Fig.  2  shows  a  sample  E(q)  with  the  fit  from  Eq.  2.  The  fit  shows  good  agreement  with  E(q)  (R2 
>0.98).  For  display  purposes,  a  one-compartment  fit  (Eq.l)  is  also  shown  to  illustrate  its  poor  agreement. 

Fig.  3  shows  bar  graphs  of  fE,  ZE,  fI?  and  Zx  fitting  results  for  each  ROI  averaged  over  all  five  SCs.  ZY  falls 

within  l-2pm,  which  is  the  expected  range  of 
axon  diameters  in  mammals8.  ZE  is  lower  than 
that  expected  for  free  water  (~20pm  for  A  = 

98.7ms).  The  ADC  calculated  from  ZE 
(-0.25x1 0'3  mm2/s)  agrees  with  literature  values 
for  fixed  spinal  cord  WM  tissue9,  which  provides 
further  evidence  that  ADCs  measured  at  low  b- 
values  (<2500  s/mm2)  primarily  reflect  diffusion 
in  ECS10.  The  average  f  was  -0.7,  which  falls 

within  the  expected  ICF  for  WM11.  An  ANOVA  analysis  indicated  no  significant  differ¬ 
ences  in  f\  among  the  WM  columns  as  previously  seen  in  mouse  SP5.  Paired  t- tests  indicated 
that  the  dorsal  column  ZE  and  Z\  are  significantly  smaller  than  those  of  the  ventral  WM  col¬ 
umn,  which  matches  previous  observations  of  smaller  MAD  and  increased  axon  density  in 
the  dorsal  compared  with  ventral  columns4,8,9. 

Conclusion 

This  work  demonstrates  the  feasibility  of  implementing  low  q-value  diffusion  MRI  on  a 
1.5T  scanner.  The  results  show  that  despite  violating  SGPA,  this  method  has  the  potential  to 
accurately  assess  regional  axonal  architecture  with  metrics  such  as  MAD  and  ICF. 
References:  1.  Beaulieu,  C,  NMR  Biomed,  15:435  (2002).  2.  Callaghan,  PT,  Principles  of  NMR  Micros¬ 
copy,  Oxford  University  Press  (1991).  3.  Cohen,  Y,  et  al.,  NMR  Biomed,  15:516  (2002).  4.  Ong,  HH, 
Neuroimage,  40:1619  (2008)..  5.  Ong,  HH,  Neuroimage,  51:1360  (2010).  6.  Price  W,  Concepts  Magn 
Res,  10:299  (1997).  7.  Kimmich  R,  NMR:  Tomography,  Diffusometry,  Relaxometry,  Springer- Verlag 
(1997).  8.  Williams,  PL,  et  al.,  Gray’s  anatomy,  Churchill  Livingstone  (1995).  9.  Schwartz,  ED,  et  al., 
AJNR,  26:390  (2005).  10.  Schwartz,  E.D.,  et  al.,  Neuroreport,  16:73  (2005).  11.  Sykova,  E,  et  al., 
Physiol  Rev,  88:1277  (2008).  Acknowledgements:  NIH  R21  EB003951 
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Fig.  3.  Bar  graphs  of  fE  (ECS  weighting),  ZE  (ECS  RMS 
displacement),  fi  (ICS  weighting),  and  Z\  (ICS  RMS  dis¬ 
placement)  fitting  results  Standard  deviation  bars  are  shown. 
Significant  p-values  (<0.05)  of  paired  t-tests  between  the 
different  ROIs  are  shown. 
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Fig.  1.  Sample  image  at  q  =  0  pm'1.  White 
bar  =  5  mm.  ROI  locations  are  shown  for  the 
dorsal,  ventral  and  lateral  WM  columns.  The 
bright  spot  above  the  center  top  SC  is  residu¬ 
al  surface  PBS. 


Fig.  2.  Sample  E{q)  for  a  lateral  WM  column 
ROI  (blue  diamonds)  with  one-compartment 
(green  line)  and  two-compartment  (red  line) 
fits. 
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Magnetic  resonance  imaging  has  previously  demonstrated  its 
potential  for  indirectly  mapping  myelin  density,  either  by  relaxo- 
metric  detection  of  myelin  water  or  magnetization  transfer.  Here, 
we  investigated  whether  myelin  can  be  detected  and  possibly 
quantified  directly.  We  identified  the  spectrum  of  myelin  in  the 
spinal  cord  in  situ  as  well  as  in  myelin  lipids  extracted  via  a  sucrose 
gradient  method,  and  investigated  its  spectral  properties.  High- 
resolution  solution  NMR  spectroscopy  showed  the  extract  compo¬ 
sition  to  be  in  agreement  with  myelin's  known  chemical  make-up. 
The  400-MHz  nH  spectrum  of  the  myelin  extract,  at  20  °C  (room 
temperature)  and  37  °C,  consists  of  a  narrow  water  resonance 
superimposed  on  a  broad  envelope  shifted  ~3.5  ppm  upfield,  sug¬ 
gestive  of  long-chain  methylene  protons.  Superimposed  on  this 
signal  are  narrow  components  resulting  from  functional  groups 
matching  the  chemical  shifts  of  the  constituents  making  up  myelin 
lipids.  The  spectrum  could  be  modeled  as  a  sum  of  super-Lorent- 
zians  with  a  T2*  distribution  covering  a  wide  range  of  values 
(0.008-26  ms).  Overall,  there  was  a  high  degree  of  similarity  be¬ 
tween  the  spectral  properties  of  extracted  myelin  lipids  and  those 
found  in  neural  tissue.  The  normalized  difference  spectrum  had 
the  hallmarks  of  membrane  proteins,  not  present  in  the  myelin 
extract.  Using  3D  radially  ramp-sampled  proton  MRI,  with  a  com¬ 
bination  of  adiabatic  inversion  and  echo  subtraction,  the  feasibility 
of  direct  myelin  imaging  in  situ  is  demonstrated.  Last,  the  inte¬ 
grated  signal  from  myelin  suspensions  is  shown,  both  spectroscop¬ 
ically  and  by  imaging,  to  scale  with  concentration,  suggesting  the 
potential  for  quantitative  determination  of  myelin  density. 

myelin  in  situ  |  myelin  NMR  spectrum  |  super-Lorentzian  fitting  |  ultrashort 
echo  time 

Myelin  is  a  critical  feature  of  nervous  system  white  matter 
(WM)  and  accounts  for  14%  of  the  wet  mass  of  WM  (1).  It 
is  a  lipid-protein  bilayer  that  extends  from  the  outer  membrane 
of  glial  cells  (i.e.,  oligodendrocytes  in  the  CNS)  and  discretely 
winds  around  individual  axonal  fibers,  leading  to  an  increase  in 
conduction  velocity  (1).  By  speeding  conduction  and  reducing 
axonal  energy  requirements,  myelin  makes  large  and  complex 
organisms  possible.  Myelin  also  contributes  to  the  mechanical 
and  functional  structure  of  the  axon.  In  addition,  some  oligo- 
dendrocytic  cells  and  precursors  can  support  action  potentials 
themselves  (2).  Deficiencies  of  myelin  lay  at  the  core  of  nu¬ 
merous  neurodegenerative  disorders,  such  as  multiple  sclerosis 
and  schizophrenia  (1).  These  deficiencies  may  result  from  develo¬ 
pmental  or  acquired  abnormalities  in  oligodendrocyte  function, 
which  also  leads  to  axonal  degeneration.  Assessment  of  myelin  may 
reveal  CNS  abnormalities  far  beyond  those  associated  with  classic 
demyelinating  diseases.  MRI  of  myelin  has  the  potential  to  char¬ 
acterize  not  only  loss  of  this  important  component  of  the  CNS  but 
also  to  reveal  axonal  and  supporting  glial  integrity  and  function. 

A  diverse  assortment  of  experimental  techniques  has  been 
applied  toward  the  goal  of  observing  and  quantifying  myelin.  The 
common  methods  rely  on  optical  microscopy  of  histologically 


stained  tissue  samples  (3).  X-ray  diffraction  (4)  and  nonlinear 
optical  techniques  (5,  6)  also  provide  insight  into  myelin  ultra¬ 
structure.  Unfortunately,  all  these  techniques  are  destructive  and 
thus  applicable  only  to  animal  studies. 

More  recently,  myelin-specific  chemical  contrast  markers  that 
selectively  bind  to  myelin  have  emerged.  Such  agents  are  cur¬ 
rently  under  development  for  both  MRI  (7)  and  positron-emis¬ 
sion  tomography  (8).  Although  these  techniques  are  potentially 
promising,  concerns  over  toxicity  may  pose  significant  hurdles  to 
their  clinical  implementation. 

So  far,  MRI  has  had  the  greatest  impact  toward  nondestruc¬ 
tive  myelin  assessment  in  both  laboratory  animals  and  humans. 
Further,  MRI  has  the  added  benefit  that  signal  contrast  origi¬ 
nates  from  endogenous  protons  and  hence  is  not  reliant  upon 
injectable  chemical  probes  nor  limited  by  contrast-related  tem¬ 
poral  delays. 

To  date,  two  indirect  MR  techniques  applicable  to  studies 
in  vivo  have  demonstrated  histologically  correlated  sensitivity 
to  myelin:  magnetization  transfer  (MT)  and  T2  relaxometry.  In 
MT,  cross-relaxation  between  myelin  protons  and  tissue  water 
is  exploited  (9).  The  signal  attenuation  resulting  from  off-reso- 
nance  saturation  (MT  ratio)  has  been  found  to  scale  with  myelin 
concentration  (10).  T2  relaxometry  yields  T2  spectra,  typically  by 
inversion  of  the  Carr-Purcell  echo  decay  using  an  inverse  Lap¬ 
lace  transformation  (11).  Spectral  components  with  T2  values 
ranging  from  10  to  50  ms  have  been  assigned  to  motionally  re¬ 
stricted  myelin  water  (12,  13)  and  have  demonstrated  strong 
correlation  with  myelin-specific  histologic  staining  (12,  14). 

Although  MT  and  T2  relaxometry  have  shown  promise,  they 
both  rely  on  indirect  detection  of  myelin  through  the  interaction 
of  water  with  myelin.  This  complex  interaction  is  affected  by 
nonmyelin  loss-related  changes,  which  can  lead  to  ambiguities  in 
data  interpretation.  For  example,  MT  is  sensitive  not  only  to 
myelin  content  but  also  to  axon  density  (15).  Therefore,  even 
though  both  techniques  may  distinguish  normal  from  abnor¬ 
mal  WM,  they  rely  on  the  invariance  of  the  myelin-water 
interaction. 

Direct  detection  of  myelin  with  MR  would  remove  some  com¬ 
plications  in  the  analysis  from  the  intermediate  effects  of  the 
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interaction  of  water  and  myelin,  and  therefore  may  provide  con¬ 
trast  specific  to  myelin  concentration.  However,  direct  detection  is 
complicated  by  motional  restriction  of  the  lipid  chains  in  the  my¬ 
elin  bilayer,  resulting  in  broad  lines  and,  consequently,  short  life¬ 
time  of  the  MR  signal. 

Perhaps  the  first  characterization  of  the  NMR  spectral  prop¬ 
erties  of  myelin  was  by  Lecar  et  al.  (16),  who  studied  anhydrous 
preparations  by  wide-line  proton  spectroscopy,  concluding  that 
under  these  conditions  the  system  is  in  a  liquid-crystalline  state. 
To  the  best  of  our  knowledge,  the  first  report  of  myelin  proton 
transverse  relaxation  was  by  Ramani  et  al.  (17).  The  authors 
performed  a  multiexponential  fit  of  spin-echo  decays  on  fixed 
human  WM  samples  from  normal  and  multiple  sclerosis  patients 
and  reported  a  T2  value  of  ~50  [is  for  myelin  protons.  They  were 
also  able  to  detect  lipid  signals  by  magic-angle  spinning  proton 
NMR  in  slices  of  normal  WM  tissue  but  not  in  multiple  scle¬ 
rosis  lesions.  Recently,  Horch  et  al.  (18)  investigated  the  T2* 
and  T2  relaxation  characteristics  of  myelin  and  reported  values 
of  ~70  [is  as  well  as  a  broader  distribution  ranging  from  50  to 
1,000  [is. 

The  transverse  relaxation  properties  of  myelin  suggest  the 
need  for  ultrashort  echo  time  (UTE)  MRI  methods,  which  entail 
collection  of  the  free-induction  decay  immediately  after  excita¬ 
tion.  Typical  implementations  include  either  3D  radial  sampling 
with  nonselective  rf  pulses  (19)  or  2D  radial  sampling  with  slice- 
selective  half  rf  pulses  and  ramp  sampling  (20)  as  previously 
applied  for  the  characterization  of  cortical  bone  matrix  and  bone 
water  (21,  22). 

UTE  MRI  has  been  used  to  image  the  short  T2*  (i.e.,  <1  ms) 
signal  from  human  brain  in  vivo  (23).  Unlike  applications  to  study 
bone,  these  implementations  include  long  T2*  suppression  meth¬ 
ods  to  attenuate  the  tissue  water  signal.  Tissue  water,  because 
of  its  rotational  mobility  and  high  concentration,  has  an  intense 
long  T2*  signal  that,  without  suppression,  overwhelms  signal  from 
short  T2*  components  (Fig.  SL4).  Although  the  images  indicated 
the  short  T2*  signal  to  be  predominantly  located  in  WM,  no  evi¬ 
dence  was  provided  to  link  this  signal  to  myelin. 

In  this  work  we  examine  the  origin  and  nature  of  the  short  T2* 
signal  of  CNS  tissue  in  freshly  excised  rat  spinal  cord  (SC)  in 
comparison  with  purified  myelin  lipid  extract  with  multinuclear 
NMR.  We  further  explore  the  potential  for  direct  detection  and 
quantification  of  myelin  by  UTE  MRI  and  discuss  the  possibili¬ 
ties  and  technical  hurdles  associated  with  translating  MRI-based 
quantification  of  myelin  to  the  clinic. 
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Fig.  1.  The  41 H  NMR  spectra  from  rat  thoracic  SC  after  D20  exchange  of  tissue 
water  (Left)  and  myelin  lipid  extract  suspended  in  D20  (Right),  showing  re¬ 
markable  similarity.  Insets:  Wide  tails  present  in  both  spectra.  Note  that  in 
addition  to  narrow  resonances,  likely  stemming  from  proteins,  the  residual 
monodeuterium  oxide  (HDO)  resonance  is  stronger  in  the  tissue  spectrum. 


Results 

High-Resolution  41 H  NMR  Spectra  of  Intact  SC  Tissue  and  Bovine  Myelin 
Extract.  Fig.  1  shows  a  comparison  of  the  NMR  spectra  col¬ 
lected  for  a  section  of  excised  rat  thoracic  SC  and  a  deuterium 
oxide  (D20)  suspension  of  myelin  lipid  extract.  The  SC  was  first 
immersed  in  99.9%  D20  for  24  h  to  exchange  tissue  water  with 
D20,  and  hence  attenuate  the  bulk  proton  signal  (Fig.  SIB).  The 
protein-free  myelin  lipid  sample  was  chemically  extracted  from 
an  intact  isolated  myelin  sample,  obtained  via  a  sucrose  gradient 
method  (SI  Materials  and  Methods)  and  quantitatively  validated 
with  proton-decoupled  13C  and  31P  NMR  spectroscopy  (Table  1 
and  Fig.  S2^4  and  B). 

Both  XH  spectra  in  Fig.  1  show  a  broad  resonance  with  a  su¬ 
perimposed  narrow  water  resonance  originating  from  residual 
HDO.  Although  the  two  broad  resonances  bear  a  high  degree  of 
similarity,  a  difference  spectrum  (i.e.,  tissue  -  extract)  (Fig.  S3) 
highlights  the  distinguishing  features.  First,  the  superimposed 
fine  structure  apparent  in  the  SC  spectrum  is  consistent  with 
mobile  proteins  as  they  might  occur  in  the  cytoplasm,  for  ex¬ 
ample.  Second,  the  difference  spectrum  contains  a  broad  reso¬ 
nance  as  it  might  be  expected  from  membrane  proteins.  The 
difference  spectrum  was  generated  such  that  the  integrated  sig¬ 
nal  area  is  27.9%  of  the  total  spectrum,  as  expected  according  to 
the  known  protein  fraction  in  myelin  (details  in  SI  Results).  Small 
errors  in  the  difference  spectrum  could  arise  because  we  are 
ignoring  cytoplasmic  proteins  in  this  fraction. 

UTE  MRI  of  Intact  SC.  Fig.  2  shows  a  set  of  images  of  freshly  ex¬ 
cised  rat  SC,  obtained  with  a  3D  radial,  ramp-sampled,  dual¬ 
echo  inversion  recovery  UTE  (de-IR-UTE)  pulse  sequence 
(Fig.  S4).  Long  T2*  tissue  water  signal  was  attenuated  via  adi¬ 
abatic  inversion  and  complex  echo  subtraction.  Adiabatic  in¬ 
version  was  used  to  significantly  reduce  the  signal  intensity  from 
tissue  water,  which  would  then  be  further  attenuated  with  echo 
subtraction.  We  empirically  selected  TI  to  achieve  the  greatest 
WM  intensity  while  minimizing  gray  matter  (GM)  intensity 
signal  in  the  complex  echo  difference  image  in  accordance  with 
the  expectation  of  low  signal  in  GM  considering  its  very  low 
myelin  content.  Images  were  collected  at  both  short  (20  [is; 
Fig.  2A)  and  long  (1,200  [is;  Fig.  2 B)  TE.  The  magnitude  of  the 
complex  difference  image  and  signal  profile  (Fig.  2  C  and  D) 

Table  1.  Lipids  of  myelin  with  abbreviations  used  in  the  text 

Molar  % 


Myelin  lipid  ID  Norton*  NMR+  %  labile  1H±at 


Cholesterol 

CHOL 

44.8 

43.1 

0.94  ±  0.04 

Galactocerebroside 

GC 

17.5 

19.6 

2.20  ±  0.39 

Galactosulfatide 

GS 

2.5 

NA§ 

0.28  ±  0.051 

Phosphatidylethanolamine 

PE 

3.4 

3.9 

0.24  ±  0.05 

PE  plasmalogen 

PEpI 

11.3 

11.8 

0.70  ±  0.13 

Phosphatidylcholine 

PC 

8.0 

7.8 

0.00  ±  0.00 

PC  plasmalogen 

PCpI 

0.3 

2.0 

0.00  ±  0.00 

Sphingomyelin 

Sph 

5.2 

5.9 

0.37  ±  0.08 

Phosphatidylinositol 

PI 

0.7 

2.0 

0.19  ±  0.03 

Phosphatidylserine 

Total 

PS 

0.2 

3.9 

0.11  ±  0.01 
5.05  ±  0.79 

Comparison  of  average  bovine  myelin  lipid  molar  ratios  reported  by  Nor- 

ton  et  al.  (25)  and  quantitative  multi-NMR  methods.  Also  shown  are  average 
percentages  of  labile  protons  pertaining  to  each  lipid  component.  PC,  phos¬ 
phatidylcholine;  PE,  phosphatidylethanolamine. 

*From  Norton  et  al.  (23). 

Present  study,  13C  and  31 P  NMR. 

^Variability  (a,  SD)  due  to  lipid  chain  length  [CH3(CH2)n;  n  =  10-25]. 

§Not  measured  owing  to  a  lack  of  an  unambiguous  resonance. 

11  Assuming  a  GS  molar  percentage  of  2.5%. 


9606  |  www.pnas.org/cgi/doi/10.1073/pnas.1115107109 


Wilhelm  et  al. 


Fig.  2.  The  3D  de-IR-UTE  images  from  rat  thoracic  SC  averaged  over  five 
central  slices.  Magnitude  images  obtained  for  (/A)  TE  =  20  |is,  ( B )  TE  =  1,200  |as, 
and  (C)  complex  difference  ( A  -  B )  (maximum-intensity  range  decreased  by  a 
factor  of  two  to  highlight  myelin  signal).  (D)  Intensity  profiles  across  the  three 
images  (delineated  as  red,  green,  and  yellow  lines  in  A,  B,  and  C,  respectively) 
to  show  relative  WM,  GM,  and  background  intensity.  The  most  intense  sig¬ 
nal,  present  in  the  short-  and  long-echo  profiles,  originates  from  residual 
surface  water.  WM  and  GM  are  indicated  in  A,  and  arrows  highlight  resid¬ 
ual  surface  water  in  B.  The  dark  boundary  observed  at  the  GMAA/M  and  WM/ 
surface  water  interfaces  in  both  echo  images  stems  from  partial  voluming  of 
adjacent  regions  with  different  T -i  resulting  in  destructive  interference. 


highlight  the  short  T2*  signal,  which  predominantly  results  from 
myelin  protons. 

MR  Signal  Dependence  on  Myelin  Concentration.  To  separate  the 
myelin  and  water  peaks  in  the  XH  spectra  of  the  myelin-D20 
suspensions,  the  spectra  were  modeled  as  a  weighted  sum  of  four 
super-Lorentzians  (SL)  for  myelin  resonances  representing  pro¬ 
tons  from  general  alkyl  chain  methylenes  (noncholesterol,  mostly 
from  fatty  acid  residues),  cholesterol  alkyl  chain  methylenes, 
terminal  methyls,  and  choline,  and  a  Lorentzian  for  the  HDO 
peak  (details  in  Materials  and  Methods).  Fig.  3A  shows  the  results 
from  fitting  of  the  proton  NMR  spectrum  of  purified  myelin 
suspended  in  D20.  The  fitting  results  were  virtually  identical  for 
all  myelin  concentrations.  Even  though  the  signal  envelope  is  very 
broad,  relatively  narrow  resonances  are  also  observed,  likely  due 
to  proton  pairs  aligned  with  an  average  orientation  at  the  magic 
angle  relative  to  the  static  field  (24). 

Relative  signal  fractions,  accounting  for  losses  during  excita¬ 
tion  and  acquisition,  along  with  associated  T2*  distributions  of 
the  four  SL  components,  were  combined  into  a  myelin  T2*  his¬ 
togram  (Fig.  3 B).  At  20  °C,  26.4%  of  the  total  signal  has  an  ef¬ 
fective  lifetime  of  <25  ps,  51.8%  of  <100  ps,  and  91.6%  of 
<1,000  ps.  At  37  °C  these  values  are  16.9%,  44.8%,  and  86.3%, 
respectively. 

Fig.  44  shows  a  series  of  fitted  myelin  signals  as  a  function  of 
decreasing  myelin  concentration.  The  NMR  signal  areas  for  the 
total  and  separate  spectral  components  (i.e.,  HDO  and  myelin) 
are  plotted  in  Fig.  4 B,  indicating  linear  scaling  with  myelin 
concentration  ( R 2  =  0.99).  We  attribute  the  positive  correlation 
of  the  water  peak  area  with  myelin  concentration  as  resulting 
from  labile  protons  from  myelin  constituents  exchanging  with 
D20  to  form  HDO.  The  calculated  average  percentage  of  labile 
protons,  for  each  of  the  10  myelin  lipid  components,  is  listed  in 
Table  1.  The  average  signal  contribution  from  the  0.1%  impurity 
of  the  D20  solvent,  calculated  as  the  y-intercept  from  the  line  of 
best  fit  (Fig.  4 B),  was  subtracted  from  all  of  the  HDO  points, 
yielding  an  estimate  of  the  labile  myelin  proton  signal  contri¬ 
bution.  The  predicted  range  of  signal  contributions  from  labile 
protons  (5.05%  ±  0.79%)  agreed  well  with  the  experimental 
HDO  peak  areas  (5.13%  ±  2.00%).  Given  the  excess  D20  used 
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Fig.  3.  'H  NMR  spectra  and  analysis  of  purified  bovine  myelin  extract  sus¬ 
pended  in  D20.  (/A)  NMR  spectrum  (black)  and  SL  fitting  showing  the 
resulting  myelin  (red)  and  HDO  (blue)  fits,  as  well  as  the  four  individual  SL 
components  of  myelin  (shaded).  The  four  fitted  SL  components  consisted  of 
(/)  SL  containing  74.3%  of  the  intensity,  centered  at  1.5  ppm,  corresponding 
to  the  general  alkyl  chain  methylene  protons  (CH2),  (/'/')  SL  containing  12.4% 
of  the  intensity,  centered  at  0.9  ppm,  corresponding  to  the  terminal  methyl 
protons  (CH3),  (Hi)  SL  containing  1 1 .1  %  of  the  intensity,  centered  at  3.2  ppm, 
corresponding  to  the  choline  methyl  protons  (Choline),  and  (/V)  SL  contain¬ 
ing  2.1%  of  the  intensity,  centered  at  1.3  ppm,  corresponding  to  the  cho¬ 
lesterol  alkyl  chain  methylene  protons  (CHOL  CH2).  (B)  J2*  histogram  of 
myelin  components  at  20  and  37  °C  derived  from  the  SL  fitting.  There  are 
small  T2*  components  that  extend  up  to  26  ms,  but  these  cannot  be  ob¬ 
served  at  the  displayed  scale.  Inset:  Myelin  extract  spectra  collected  at  the 
two  temperatures. 


in  the  suspensions,  it  is  reasonable  to  assume  that  all  labile 
myelin  protons  had  exchanged  with  deuterium. 

Fig.  5  Inset  shows  the  complex  difference  2D  projection  de-IR- 
UTE  image  for  a  series  of  myelin  suspensions  of  increasing 
concentration.  Region  of  interest  (ROI)  average  intensities  from 
each  of  the  myelin  samples  in  the  image  are  plotted  in  Fig.  5  and, 
analogous  to  the  spectral  data,  are  linearly  correlated  with  my¬ 
elin  concentration  ( R 2  =  0.98). 

On  the  basis  of  the  relaxation  characteristics  of  myelin  (notably 
the  lifetime  of  the  various  T2*  components),  Bloch  equation 
simulations  (details  in  SI  Materials  and  Methods)  accounting  for 
losses  sustained  during  the  rf  pulse  sampling  suggest  that,  under 
the  spectral  recording  conditions,  81.7%  of  the  total  spectral  sig¬ 
nal  is  recovered,  where  the  shorter  T2*  components  account  for 
the  majority  of  the  signal  lost.  The  depletion  of  the  imaging  signal 
is  more  severe  because  it  entails  coherence  losses  during  both 
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Fig.  4.  04)  Fitted  myelin  peaks  for  myelin  suspensions  of  various  concen¬ 
trations.  ( B )  Linear  correlation  plot  of  MR  signal  as  a  function  of  myelin  con¬ 
centration  for  the  total  (triangle),  myelin  (circle),  and  HDO  (square)  signal 
components  with  calculated  R2  values.  All  two-tailed  P  values  <0.01.  Also 
shown  is  the  calculated  signal  fraction  of  labile  myelin  protons  (green  shaded). 

the  adiabatic  inversion  and  excitation  pulses,  as  well  as  during  the 
period  preceding  collection  of  the  low-spatial  frequencies  (TE) 
and  partial  saturation.  We  estimate  that  with  the  imaging  param¬ 
eters  used  in  our  experiments,  38.0%  of  the  total  available  myelin 
proton  signal  could  be  recovered. 

Discussion 

In  this  work  we  explored  the  feasibility  of  direct  imaging  and 
quantification  of  myelin  by  magnetic  resonance  as  an  alternative 
to  indirect  detection  techniques  such  as  MT  or  T2  relaxometry. 
Direct  detection  of  myelin  would  remove  any  ambiguities  in 
analysis  and  provide  a  contrast  specific  to  myelin  concentration. 
To  explore  the  feasibility  of  direct  myelin  imaging,  we  first  in¬ 
vestigated  the  origins  of  the  short  T2*  signal  in  intact  SC  and 
myelin  extracts.  We  then  presented  preliminary  UTE  with  long 
T2*  suppression  images  of  myelin  in  WM. 

The  purity  of  our  myelin  extract  was  validated  by  9.4  T  high- 
resolution  proton-decoupled  31P  (Fig.  S24)  and  13C  (Fig.  S2 B) 
NMR,  showing  good  agreement  with  previously  determined  mo¬ 
lar  ratios  (25)  (Table  1).  The  slight  variations  are  reasonable  given 
the  differences  in  the  anatomical  origin  of  the  sample  (brain  vs. 
SC)  (26),  inherent  variability  of  the  measurement  (27),  and  the 
developmentally  immature  nature  of  the  tissue  examined  (26).  Our 
13C  spectra  also  agreed  well  with  Husted  et  al.’s  reported  proton- 
decoupled  magic-angle  spinning  (MAS)  13C  NMR  spectra  of 
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Fig.  5.  UTE  imaging  and  analysis  of  purified  bovine  myelin  extract  sus¬ 
pended  in  D20.  Inset:  The  2D  projection  de-IR-UTE  image  of  six  5-mm  NMR 
tubes  filled  with  myelin/D20  suspensions  of  varying  concentrations  (mg/mL). 
Mn  indicates  the  tube  containing  manganese-doped  water.  Plot  of  the 
normalized  de-IR-UTE  mean  ROI  intensities  vs.  myelin  concentration  (P  = 
0.0002).  Arrows  indicate  areas  of  signal  from  paper  towel  used  to  stabilize 
the  tubes. 

isolated  human  myelin  (28).  In  that  study  a  synthetic  13C  spectrum 
of  myelin  was  generated  as  a  weighted  sum  of  individual  lipid 
spectra,  based  upon  the  known  ratios  of  the  individual  lipid  com¬ 
ponents.  The  synthetic  spectrum  was  found  to  be  in  good  agree¬ 
ment  with  the  observed  MAS  spectrum,  which  in  turn  closely 
matches  our  solution  spectrum  (Fig.  S2 B),  except  for  the  absence 
of  signals  assigned  to  protein  resonances  in  the  MAS  spectrum. 

Any  conclusions  from  extracted  or  synthetically  produced  my¬ 
elin  beg  the  question  whether  the  resulting  material  retains  the 
biophysical  properties  of  native  myelin.  It  is  known  that  in  aque¬ 
ous  suspensions  myelin  lipids  spontaneously  self-assemble  into 
bilayer  structures  (i.e.,  vesicles).  In  this  enthalpy-driven  process 
(29)  the  myelin  lipid  suspension  regains  a  structural  order  remi¬ 
niscent  of  physiological  myelin  minus  the  protein  component. 

The  spectrum  of  the  isolated,  reconstituted  myelin  exhibited  a 
very  broad  line  with  relatively  narrow  components  centered  ~3.5 
ppm  upfield  from  water,  consistent  with  methylene  protons  of 
alkyl  chains,  the  main  constituent  of  myelin.  Further,  this  spec¬ 
trum  bears  a  high  degree  of  similarity  with  that  of  neural  tissue, 
which  supports  the  notion  that,  upon  aqueous  suspension  of  the 
extract,  a  bilayer  structure  analogous  to  that  for  native  myelin  is 
reconstituted.  The  SC  spectrum  is  expected  to  contain  nonmyelin 
contributions  such  as  from  proteins.  The  difference  spectrum 
(Fig.  S3)  suggests  the  presence  of  motionally  restricted  mem¬ 
brane  proteins,  along  with  a  small  fraction  of  mobile  (perhaps 
cytoplasmic)  proteins  yielding  narrow  resonances. 

The  myelin  resonance  is  consistent  with  the  SF  lineshape  of 
a  dipolar-broadened  liquid-crystalline  lipid  system  (24,  30)  as 
suggested  by  several  lineshape  properties  (31),  including  a  small 
second  moment  (M2  =  1.06  x  10g  s-2),  a  large  fourth  moment 
(M4)  to  M2  ratio  [M4/(M2)2  =  6.57],  and  a  full  width  at  1/10 
maximum  (Avi/i0)  to  full  width  at  1/2  maximum  (Aiq/2)  ratio  >3 
(Avi/i0:  Av4/2  =  7).  Following  previous  work  (32),  we  modeled 
the  myelin  lipid  extract  spectrum  as  a  sum  of  SFs  representing 
protons  from  general  alkyl  chain  methylenes,  cholesterol  alkyl 
chain  methylenes,  terminal  methyls,  and  choline.  Although  the 
T2*  distribution  showed  a  wide  range  of  values  (0.008-26  ms), 
the  distribution  was  dominated  by  that  of  the  alkyl  chain  meth¬ 
ylenes  because  they  make  up  ~75%  of  the  four  myelin  lipid 
proton  moieties  under  consideration.  One  notes  a  bimodal  dis¬ 
tribution  (Fig.  S5)  with  the  first  peak  significantly  attenuated  as 
expected  for  an  SF  lineshape  determined  by  the  angular  de¬ 
pendence  of  T2*  (details  in  SI  Results). 

In  contrast  to  the  SF  fitting  performed  in  this  work,  Horch 
et  al.  performed  multiexponential  fitting  of  the  time-domain 
signal  decay  (18).  Although  the  SF  lineshape  theoretically 
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cannot  be  described  with  a  basis  set  of  exponential  functions,  the 
authors  suggested  that  the  resulting  errors  would  be  small  for  the 
case  of  myelin.  Horch  et  al.  analyzed  free-induction  decay  signals 
of  myelinated  mammalian  and  amphibian  nerves,  as  well  as 
synthetic  myelin  at  4.7  T,  yielding  histograms  of  relaxation  times. 
The  authors  detected  significant  components  with  T2*  values  of 
20  and  70  ps  in  frog  sciatic  nerve,  which  they  conjectured  to  arise 
from  protein  and  methylene  protons  of  myelin,  respectively.  In 
contrast,  because  our  myelin  extract  was  free  of  protein,  the 
present  data  alternatively  suggest  the  short  T2*  components 
(<25  ps)  to  arise  from  myelin  lipids.  This  is  not  to  imply  that 
membrane  proteins  cannot  contribute  a  short  T2*  component,  as 
suggested  by  the  broad  resonance  in  the  difference  spectrum 
(Fig.  S4)  and  in  other  reports  (24,  33). 

In  bovine  myelin  suspended  in  D20,  we  found  the  integrated 
spectral  areas  to  scale  linearly  with  myelin  concentration  (Fig.  4), 
as  did  ROI  intensities  of  the  2D  projection  de-IR-UTE  images 
(Fig.  5),  thus  suggesting  that  quantitative  myelin  imaging  may  be 
feasible.  Direct  3D  de-IR-UTE  imaging  of  a  rat  SC  in  situ  at  400 
MHz  highlights  the  potential  of  such  an  approach,  as  demon¬ 
strated  with  images  showing  signal  from  the  WM  regions  only. 
Absolute  quantification  would  require  a  reference  sample,  ide¬ 
ally  with  relaxation  and  density  properties  matching  myelin.  The 
reference  sample  should  also  be  chemically  stable. 

Previously,  Waldman  et  al.  (23)  obtained  images  of  the  human 
brain  using  a  slice-selective  UTE  along  with  soft-tissue  suppres¬ 
sion,  essentially  showing  intense  signal  from  WM  regions  of 
the  brain,  which  they  attributed  to  short  T2*  components.  Under 
these  conditions  [i.e.,  selective  half -sine  pulses  of  400-  to  600-ps 
duration  (34)],  all  but  the  longest  T2*  components  of  the  myelin 
protons  would  elude  detection. 

Our  results  indicate  approximately  20%  of  the  protons  in 
myelin  lipids  to  have  an  effective  T2*  less  than  25  ps.  Even  under 
the  more  favorable  conditions  of  our  imaging  experiments,  the 
majority  of  these  short  T2*  components  remains  undetectable. 
The  very  short  lifetime  of  the  signal  has  potentially  adverse 
effects  on  the  point-spread  function  (PSF)  manifesting  as  blur¬ 
ring.  A  simple  estimation  based  on  the  FWHM  of  the  spectra 
(Fig.  3)  predicts  an  intrinsic  resolution  (defined  as  the  minimum 
achievable  pixel  size)  on  the  order  of  100  pm,  which  is  approx¬ 
imately  one  pixel  with  our  current  imaging  parameters  (more 
details  in  SI  Materials  and  Methods).  It  should  be  noted  that  al¬ 
though  the  blurring  from  the  shortest  T2*  component  would  be 
greater,  its  impact  on  the  PSF  is  negligible  because  this  signal 
fraction  remains  virtually  undetected. 

A  limitation  that  needs  to  be  noted  for  this  method  is  that  it 
detects  myelin  solely  on  the  basis  of  its  T2*  properties.  Thus, 
errors  in  long  T2*  suppression  may  lead  to  signal  misclassified  as 
short  T2*  and  hence  falsely  identified  as  myelin.  Such  errors 
could  be  accounted  for  and  perhaps  mitigated  by  tailoring  a 
reference  sample  so  as  to  contain  water  of  comparable  concen¬ 
tration  and  relaxation  times  to  those  of  neural  tissue.  Last,  there 
are  other  possible  nonmyelin  short  T2*  sources  that  could  con¬ 
tribute  to  the  UTE  image  intensity,  including  glial  cell  mem¬ 
branes,  calcifications,  tissue  scars,  vasculature,  and  hemorrhage 
(23),  that  would  be  indistinguishable  from  myelin. 

Another  potential  problem  could  arise  from  saturation  of  the 
myelin  signal  via  cross-relaxation  (35).  Even  though  adiabatic  in¬ 
version  of  tissue  water  has  minimal  direct  effect  on  the  myelin 
lipid  proton  signal,  transfer  of  magnetization  from  the  water  to  the 
myelin  proton  pool  could  occur  via  chemical  exchange  or  dipolar 
coupling.  This  mechanism  would  result  in  a  reduction  in  the  de¬ 
tected  myelin  signal,  an  effect  that  requires  further  scrutiny. 

The  potential  for  translation  of  the  method  to  the  clinic  will  be 
challenging.  Nevertheless,  it  is  encouraging  to  note  that  with 
dedicated  hardware  rf  pulses  of  20  ps  or  less  have  already  been 
shown  to  be  feasible  on  clinical  equipment  at  3  T  field  strength, 
as  in  recent  work  by  Wu  et  al.,  who  imaged  the  collagen  matrix 


of  cortical  bone  (36).  Further,  at  37  °C  the  measured  T2*  values 
increase  by  ca.  30%  (Fig.  3 B),  raising  the  minimum  T2*  value 
from  8  to  10.5  ps.  At  body  temperature  the  scan  parameters  used 
in  our  current  experiments  at  9.4  T  (except  for  a  20-ps  hard  pulse 
of  7.6°  flip  angle  to  match  the  peak  Bx  amplitudes  of  clinical  head 
coils)  predict  4.9%  of  the  total  myelin  proton  signal  [i.e.,  0.7%  of 
the  total  proton  signal  given  that  myelin  accounts  for  14%  of  WM 
(1)]  to  be  recoverable  on  a  3  T  clinical  MRI  system  (for  calcu¬ 
lation  details  see  SI  Materials  and  Methods). 

Given  that  tissue  proton  concentration  is  ~100  M,  the  con¬ 
centration  of  detectable  myelin  protons  is  approximately  700  mM. 
In  comparison  with  proton  MR  spectroscopic  imaging  (MRSI)  of 
brain  metabolites,  where  the  metabolite  concentrations  are  on 
the  order  of  10  mM,  detectable  myelin  proton  concentration, 
and  hence  intrinsic  SNR,  should  be  a  factor  of  approximately  70 
greater  than  that  of  typical  metabolites.  However,  this  gain  in 
SNR  compared  with  MRSI  is  mitigated  by  the  reduced  sampling 
time  imposed  by  the  much  shorter  T2*  of  the  protons  in  myelin 
compared  with  those  in  brain  metabolites.  We  estimate  reduced 
overall  sampling  time  to  result  in  a  loss  on  the  order  of  a  factor 
of  10.  Given  a  reported  resolution  for  MRSI  of  5-10  mm  (37), 
we  project  the  resolution  achievable  with  our  method  to  be 
roughly  on  the  order  of  2.5-5  mm  with  T2*-induced  PSF  blurring 
not  exceeding  0.6  mm  (SI  Results). 

Conclusions 

We  have  characterized  the  spectral  properties  of  the  myelin 
proton  signal  in  situ,  as  well  as  in  reconstituted  suspensions  of 
myelin  lipid  extract.  Our  results  show  that  the  short  T2*  com¬ 
ponent  of  WM  originates  primarily  from  myelin  lipid  protons 
and  further  that  direct  imaging  of  these  protons  is  possible  even 
though  the  shortest  components  are  not  detectable.  Last,  our 
analysis  suggests  that  translation  from  the  laboratory  to  clinical 
MRI  will  be  challenging. 


All  MR  spectroscopy  and  imaging  experiments  were  performed  on  a  9.4  T 
vertical-bore  spectrometer/microimaging  system  (Bruker  DMX  400)  with 
Micro2.5  gradients  (100  G/cm  maximum  strength)  and  BAFPA40  amplifiers. 

Neural  Tissue  Preparation.  SC  samples  were  harvested  from  healthy  adult 
Sprague-Dawley  rats  (Charles  River  Laboratories)  and  bovine  spinal  columns 
(Bierig  Brothers  Veal  and  Lamb  Products).  The  rats  were  killed  by  carbon 
dioxide  asphyxiation  in  accordance  with  an  Institutional  Animal  Care  and  Use 
Committee-approved  protocol.  After  killing,  rat  spinal  columns  were  re¬ 
moved,  and  the  SC  was  dissected  out. 

NMR  Spectroscopy.  High-resolution  proton-decoupled  13C  NMR  [Sweep  width 
(SW)  =  24  kHz,  number  of  scans  (NS)  =  36,768,  number  of  real  and  imaginary 
data  points  (TD)  =  65,536,  repetition  time  (TR)  =  1.36  s,  a  =  30°]  and  proton- 
decoupled  31P  (SW  =  3.23  kHz,  NS  =  8,536,  TD  =  8,192,  TR  =  1.27  s,  a  =  30°) 
spectra  were  collected  for  samples  of  purified  bovine  myelin  extract,  dis¬ 
solved  in  a  (5:4:2)  ternary  mixture  of  deuterated  chloroform  (99.8  atom  %  D; 
Acros  Organics),  methanol  (99.8  atom  %  D;  Acros  Organics),  and  0.2  M  EDTA/ 
water  (99.9  atom  %  D;  Sigma-Aldrich). 

All  NMR  spectra  were  collected  with  the  following  parameters:  SW  = 
100  kHz,  NS  =  256,  TD  =  262,144,  TR  =  3.6  s,  a  =  90°,  pulse  duration  =  9.6  ps. 
Freshly  excised  SC  sections  (<2-h  postmortem  interval)  were  immersed  in  a 
perfluorinated  oil  (Fomblin-Y;  Sigma-Aldrich)  before  experiments. 

SL  Fitting  of  Proton  Spectrum.  As  described  by  Wennerstrom  (30),  only  partial 
averaging  of  dipolar  coupling  via  translational  and  rotational  diffusion 
occurs,  resulting  in  a  dipolar-broadened  liquid-crystalline  lipid  system  with 
an  SL  lineshape  that  can  be  written  as: 

«*>-/l55J§b I'lsSiM1* 

where  oo0  is  the  chemical  shift,  0  is  the  angle  of  the  lipid  bilayer  surface 
normal  with  respect  to  B0,  and  f  (cd-cd0)  is  any  highly  peaked  lineshape  such  as 
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a  Gaussian  or  Lorentzian.  Assuming  0  is  uniformly  sampled  and  setting  f 
(< ro-a)o )  to  be  a  Lorentzian,  it  can  be  seen  from  Eq.  1  that  an  SL  is  composed 
of  a  series  of  scaled  Lorentzians.  From  the  widths  and  intensities  of  these 
Lorentzians,  the  T2*  distribution  of  a  single  SL  can  be  calculated.  Protons  at 
different  chemical  shifts  (e.g.,  alkyl  chain  methylenes,  terminal  methyls,  and 
choline)  are  each  expected  to  give  rise  to  SL  lineshapes  (32). 

Spectral  fitting  was  performed  in  Matlab  (Mathworks).  Four  SLs  were  used 
to  represent  general  alkyl  chain  methylenes  (noncholesterol,  mostly  from 
fatty  acid  residues),  cholesterol  alkyl  chain  methylenes,  terminal  methyls,  and 
choline,  whereas  a  single  Lorentzian  was  used  to  model  residual  HDO.  Be¬ 
cause  cholesterol  alkyl  chain  methylenes  sit  deep  within  the  lipid  bilayer,  it  is 
reasonable  to  expect  them  to  be  more  mobile  than  the  general  alkyl  chain 
methylenes,  therefore  resulting  in  a  narrower  SL.  The  chemical  shifts  of  each 
SL  were  set  to  the  known  isotropic  shift  of  the  various  moieties,  and  the  width 
and  relative  intensities  were  free  parameters.  The  R2  of  the  fit  was  greater 
than  0.99. 

UTE  MR  Imaging.  The  3D  de-IR-UTE  imaging  (Fig.  S4):  SW  =  200  kHz,  TE  = 
20/1,200  |is,  Tl  =  500  ms,  TR  =  1  s,  field  of  view  =  15  mm,  matrix  size  =  128  x 
128  x  128,  number  of  views  =  5,342,  pulse  duration  =  20  ps.  The  sequence 
was  based  on  that  used  by  Anumula  et  al.  (38).  Tl  was  determined  empiri¬ 
cally  as  the  duration  yielding  optimal  GM  suppression  (because  GM  is  ex¬ 
pected  to  have  negligible  myelin  concentration)  in  a  complex  difference 
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Introduction 

Deficiencies  of  myelin,  a  lipid  bilayer  sheath  critical  for  normal  function  of  white  matter  (WM),  lay  at  the  core  of  numerous  neurodegenerative 
disorders  such  as  multiple  sclerosis  and  schizophrenia  (1).  At  present,  there  are  few  alternatives  to  destructive  histologic  methods  to  directly  assess 
myelin.  The  short  T2  of  myelin  protons  make  ultra-short  echo  time  (UTE)  MRI  a  potential  imaging  modality  to  directly  detect  myelin  (2).  In 
contrast,  indirect  MRI  methods  such  as  magnetization  transfer  and  T2  relaxometry  are  based  on  complex  interactions  between  water  and  myelin, 
which  can  lead  to  ambiguities  in  data  analysis.  Characterizing  the  T2*  distribution  of  myelin  is  key  to  developing  optimal  UTE  methods  for  myelin 
imaging.  Previous  attempts  have  used  multi-exponential  fitting  of  the  FID  (3),  which  is  not  only  an  ill-posed  problem  (4),  but  also  theoretically 
incorrect.  Myelin  is  a  liquid  crystalline  lipid  system  that  is  described  by  a  sum  of  super-Lorentzians  (SL)  rather  than  a  multi-Lorentzian  lineshape  (5, 
6).  Here,  we  use  this  SL  framework  to  calculate  T2  distributions  by  fitting  !H  NMR  spectra  of  myelin  lipid  extract  and  intact  rat  spinal  cord  (SC). 
Theory  and  Methods 

According  to  Wennerstrom  (6),  due  to  averaging  effects  from  translational  and  rotational  diffusion,  for  a  given  orientation  of  a  lipid  bilayer,  the 
lineshape  can  be  expressed  as:  L(  5-50,  0)  =  |3cos20  -If1/  [(5-50)/  |3cos20  -1|  ]  (Eq.l),  where  5  is  the  chemical  shift  centered  at  80,  0  is  the  angle  of  the 
lipid  bilayer  surface  normal  with  B0,  and  f(x)  is  any  highly  peaked  lineshape  such  as  a  Gaussian  or  Lorentzian.  The  SL  lineshape,  Lsl(8-80),  results 
from  a  uniform  sampling  of  0  from  0  to  7i/2:  Lsl(8-80)  =  /L(8-80,  0)  sin0  d0  (Eq.2).  By  assuming/(v)  to  be  a  Lorentzian,  it  can  be  seen  from  Eqs.  1 
and  2  that  a  SL  is  composed  of  a  series  of  scaled  Lorentzians.  From  the  widths  and  intensities  of  these  Lorentzians,  the  T2  distribution  of  a  single  SL 
can  be  calculated.  Multiple  SLs  have  been  used  to  fit  NMR  spectra  of  model  membrane  systems  in  which  the  SLs  arise  from  protons  at  different 
chemical  shifts,  e.g.  alkyl  methylenes,  terminal  methyls,  and  choline  (7).  Therefore,  it  is  possible  in  theory  to  perform  a  multi-SL  fit  of  a  lH  NMR 
spectrum  of  myelin  and  calculate  a  T2*  distribution. 

Rat  and  bovine  SC  samples  were  harvested  from  Sprague-Dawley  rats  (Charles  River  Labs)  and  a  local  butcher.  Myelin  lipids  were  extracted 
from  bovine  SC  tissue  with  a  sucrose  gradient  method  (8),  dissolution  in  a  ternary  mixture  (chloroform/methanol/water),  and  lyophilization.  Previous 
work  has  shown  that  this  protocol  extracts  myelin  lipids  with  little  to  no  protein  (2).  Dehydrated  myelin  lipid  extract  was  then  re-suspended  in  99.9% 
D20  (Sigma- Aldrich)  to  regenerate  a  bilayer  structure.  lH  NMR  spectra  at  9.4T  (DMX-400,  Bruker  Instruments)  were  obtained  for  a  freshly  excised 
rat  thoracic  SC  immersed  in  Fomblin  (Sigma- Aldrich),  as  well 
as  the  myelin  lipid  extract.  Rat  SC  was  immersed  in  D20  for 
24  hrs  prior  to  experiments  to  reduce  the  dominant  tissue 
water  peak. 

Spectral  fitting  was  performed  in  Matlab  (Math works). 

Four  SLs  were  used  to  represent  general  alkyl  chain 
methylenes,  cholesterol  alkyl  chain  methylenes  (as  they  have 
shorter  chain  lengths),  terminal  methyls,  and  choline,  while  a 
single  Lorentzian  was  used  to  model  residual  HDO.  The 
chemical  shifts  of  each  SL  were  set  to  the  known  isotropic 
shift  of  the  various  moieties,  and  the  width  and  relative 
intensities  were  free  parameters. 

Results  and  Discussion 

Fig.  1  shows  myelin  lipid  extract  and  rat  SC  !H  NMR  spectra  with  the  results  of  the  four-component  SL  fits.  Both  spectra  are  comprised  of  a 
narrow  peak  from  residual  HDO  and  a  broad  resonance  (linewidth  ca.  1700Hz)  from  myelin.  A  large  4th  to  2nd  moment  ratio  (ca.  6.6),  suggests  that 
this  broad  resonance  has  a  SL  lineshape.  The  SC  l¥L  NMR  spectrum  had  additional  minor  peaks  from  intracellular  proteins  and  other  non-myelin 
protons.  Despite  the  additional  complexity  of  the  SC  spectrum,  both  SL  fits  agree  well  with  the  !H  NMR  spectra  (R2  >  0.99).  Fig.  2  shows  the  relative 
fractions  (theoretical  and  fitted)  of  the  four  SL  myelin  components  in  myelin  extract  and  intact  rat  SC.  As  expected,  the  signal  is  dominated  (>70%) 

by  alkyl  methylene  protons.  Deviations  from  theoretical  values  may 
result  from  inaccuracies  of  the  SL  framework  to  describe  non-chain 
alkyl  protons,  e.g.  choline  and  terminal  methyls. 

Fig.  3  shows  the  T2  distributions  derived  from  the  SL  fits  for  the 
myelin  lipid  extract  and  SC.  The  T2*  distributions  are  highly  skewed 
with  a  wide  range  (lOjas  to  10ms).  Despite  this  range,  roughly  50% 
(80%)  of  the  signal  has  a  T2*  less  than  20ps  (lOOps).  This  result 
highlights  the  difficulty  of  direct  myelin  imaging  even  with  UTE  MRI. 
Further  investigation  is  needed  to  study  the  system  at  body  temperature 
as  increased  molecular  motion  is  likely  to  result  in  longer  effective  T2  s. 
Conclusion 

This  work  uses  a  SL  framework  to  characterize  the  T2  distribution 
of  myelin,  which  would  provide  guidance  toward  developing  UTE 
methods  for  myelin  imaging.  The  results  indicate  that  at  ambient 
temperature  ~50%  of  the  myelin  lipid  proton  signal  has  T2*  <  20  jis. 
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Fig  2.  Relative  fraction  fitting  results 
of  the  four  SL  components  with 
expected  theoretical  fractions.  Error 
bars  account  for  variation  in  alkyl 
chain  length. 
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Fig  3.  Calculated  T2*  distributions  for 
myelin  lipid  extract  (red)  and  rat  SC 
(blue).  Signal  fractions  with  T2*  <20  and 
<100  ps  are  reported. 


Fig  1.  lR  NMR  spectra  at  ~20°C  (black)  and  SL  based  fit  (red)  for  a)  myelin  extract 
(R2=0.999)  and  b)  rat  SC  (R2=0.997).  The  four  SL  and  Lorentzian  peaks  of  the  fit  are  shown. 
Full  spectral  width  is  shown  in  inset  with  HDO  peak  truncated. 
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